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Abstract: 

Recently, new thermodynamic inequalities have been obtained, which set bounds on the 
quadratic fluctuations of intensive observables of statistical mechanical systems in terms of 
the Bogoliubov - Duhamel inner product and some thermal average values. It was shown 



a* 



that several well-known inequalities in equilibrium statistical mechanics emerge as special 
cases of these results. On the basis of the spectral representation, lower and upper bounds 
on the one-sided fidelity susceptibility were derived in analogous terms. Here, these results 
are reviewed and presented in a unified manner. In addition, the spectral representation of 

\0 ■ the symmetric two-sided fidelity susceptibility is derived, and it is shown to coincide with 

the one-sided case. Therefore, both definitions imply the same lower and upper bounds on 

■^ ' the fidelity susceptibility. 
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I. INTRODICTION 

A number of important results on the role of critical fluctuations in systems with broken 
symmetries have been obtained in the past by using famous inequalities due to Bogoliubov, 
Mermin- Wagner, Griffits, among others, see, e.g. [lH4(. A remarkable property of this 
approach is that the results obtained are exact and cannot be inferred from any perturbation 
theory. The benefits of having exact statements about systems with strongly interacting 
particles are difficult to overestimate. 

The article is structured as follows. We begin by introducing notations and basic defi- 
nitions of the models under consideration. In Section II, we review the recently obtained 
inequalities, in which the central role is played by the Bogoliubov - Duhamel inner product. 
In Section III we give the definitions of fidelity and its second derivatives: the one-sided and 



symmetric two-sided fidelity susceptibilities. The new results concerning the derivation of a 
spectral representation for the symmetric two-sided fidelity susceptibilities are presented in 
Section IV. Since it turns out that the two spectral representations are the same, we note 
that both definitions imply the same lower and upper bounds, see Section V. The paper 
closes with some concluding remarks in Section VI. 

We shall consider families of Hamiltonians depending on auxiliary external fields v and 
u* conjugate to the intensive operators A^ and A, respectively, 

H{v)=H-\A\(uAi + v*A). (1) 

The quantum statistical mechanical models are defined initially in a finite region A of the 
d-dimensional Euclidean space R d or the d- dimensional integer lattice Z d . By |A| we denote 
the volume of A in the former case, or the number of lattice sites in the latter case. 

The Hamiltonian H-[y) is defined as a self-adjoint operator in a separable Hilbert space 
H, and the corresponding free energy density f\\H{v)\ is assumed to exist. For the sake 
of simplicity, we do not explicitly distinguish between a Hamiltonian H.[y), describing a 
system with fixed number of particles N in A, and the statistical operator *H(is) — jmN" in the 
grand canonical ensemble, where \x is the chemical potential and j\f is the particle number 
operator. The density of the corresponding thermodynamic potential is given by 

U[H{v)\ = -{p\h\)- 1 \nZ t [U{ V )l (2) 

where Z[H(is)] := Tre~^ H ^ v ' is the partition function. Average values in the Gibbs ensemble 
with the Hamiltonian "H(v) are defined as 

(■■■} n{u) = Tr(e-^...)/Z A [H(v)]. (3) 

In the theory of phase transitions a key role play the following quantities: the quadratic 
fluctuations 

(6A6A% (V) , SA = A-(A) n{vh (4) 

and the second derivative of the free energy density with respect to external fields is, is*, i.e., 
the isothermal susceptibility 

d*f A [H(is)] 

XA{U) = ~ dis+dis • (5) 

In what follows we will focus on the problem of relating (j3j) and fl5]) to other quantities 
that appear as indicators of phase transitions. Much work have been done on this issue in 
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the context of quantum information theory (see, e.g., p, |6[ and references therein). As a 
result, the knowledge obtained so far shows that techniques borrowed from the quantum 
information theory are likely to be useful in studying the phenomenon of phase transitions. 
First of all we make use of the important relationship between (JH]) and the Bogoliubov - 
Duhamel inner product, defined as 

(A; B) n = {Z[H})- 1 f dr Tr [e^^^V^S] , (6) 

Jo 

namely, 

'WMU^M). (7) 



P\A\ dv*dv 

This object has been introduced and discussed by a number of authors (see [2J, |7| and refer- 
ences therein), which noted that it is useful in finding quantum generalizations of different 
inequalities in statistical mechanics of classical systems. Let us note that if the operator 
A commutes with the Hamiltonian % (commutative case) the expressions (jl]) and (J7|) coin- 
cide. The genuine difference between (J3J) and (J7J) in the noncommutative case constitutes 
the novelty of the inequalities obtained below. 

In the remainder we consider given observables A,B,C,... and fixed Hamiltonian % 
pertaining to a quantum system in a finite region A. Whenever no confusion arises, for 
brevity of notation we will omit the subscripts "H and A, as well as the argument of the 
partition function Z. 

II. BOUNDS ON THE QUADRATIC FLUCTUATIONS 

We assume that the Hamiltonian 7i of the system in a finite region of space A is a self- 
adjoint, trace-class operator which generates the Gibbs semigroup {exp(— /3"Ha)}/3>o- Let the 
operators A, B, . . . , belong to the algebra of bounded observables, for which the Bogoliubov 
- Duhamel inner product (A; B) is well defined. 

Further we assume that the Hamiltonian "H has a discrete, non-degenerate spectrum only, 
{E n , n — 1, 2, 3 . . . } and denote by \n) the corresponding eigenfunctions, i.e., W\n) = E n \n), 
n = 1,2,3.... By A mn = (m\A\n) we denote the corresponding matrix element of an 
operator A. Then, the spectral representation of the Bogoliubov - Duhamel inner ([6]) can 



be written as 



p—f3E m p—fiE n 



(A; B) H = (Z A N)- 1 X; , ^ mB B(E n -E m ) + t^^)" 1 E e ~^An, ^ 

m,n *■ ' n 

where the prime in the double sum means that the term with n = m is excluded. 
Our aim is to majorize the quadratic fluctuations 

(6A^A) = (A^A) H -\(A) H \ 2 (9) 

by terms proportional to some power of the inner product 



(8A;5A) = (A;A)-\(A)\ 
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g -pE m „—/3E n 



^E'lAnnl * f} {En _ Em) +{Z h [U])-^- PEn \Ann\ 2 -\{A)\\ (10) 



rnji 



Fot the symmetrized form of (Q we obtain 

l -^A + AA^)-{A-A) 



, , ,1 p— 0E m p —/3E n 

^ £'MM 2 i(e-**. + e"^) - e J , (11) 



m,,n 



/3(E n — E n 



which, by using the identity 



-PE m + Q -pE n = (e -^„ _ e - P E m) CQth P(Em E n ) (i2) 



can be expressed as: 



I p — ftEm p — /3E„ 

-{A A + Atf) - (A; A) = Z- 1 V '| A mn \ 2 — — — - (X mn cothX mn - 1), (13) 

2 tZ P\^n ~ ^" 



J m I 



where X mn = (3(E m - E n )/2. 

Different choices of the upper bound on the right-hand side of ([TBI generate different 
inequalities. Thus, the inequality of Brooks Harris [8] , 

(A; A) < 1 -(AA+ + A+A) < (A;A) + ^([[A\n},A}}. (14) 

is obtained by setting 1 < x cothx < 1 + x 2 /3. 

On the other hand, if one uses another elementary inequality, 1 < xcothx < 1 + |x|, and 
subsequently applies the Holder inequality, one obtains the result due to Ginibre [9j: 

(A; A) < \{AA + + A + A) < (A; A) + 1 -{{A- A)(3([[A + ,n], A])}*. (15) 



A different choice of the parameters in the Holder inequality, followed by the implemen- 
tation of the upper bound 

I e -m _ e -? Em \ < | e - m + e -p Em ^ (16) 

generates a symmetric version of the inequality due to Bogoliubov Jr. [10]: 

l -(AA + + A + A) < (A;A) + ±{(A;A)l3} 2 / 3 {({A\H]iH,A] + {H,A}{A + ,H})} 1 / 3 , (17) 

To derive generalizations of the known inequalities involving the Bogoliubov - Duhamel 
inner product, we define a set of new functionals in terms of their spectral representations: 



\2n-l 

Wl 



F 2n (J; J) = Z- 1 J2 \Jmi\ 2 \e~ fiEl - e-P E ™\(P\E m - E l 

ml 

= (3 2n (R n , R n ) = f3 2n - 1 ([RtR n -i - R n -iR n ]), n = 0,1,2,3,. . ., (18) 

where, by definition, i?_i = Xj^ is a solution of the operator equation J = [Xjn,'H], and 
R = R (J) = J, R 1 = R 1 (J) = [H,J], ..., R n = R n (J) = [H,R n ^(J)]. (19) 



The observables Rk, k = 0, 1, 2, . . . , were introduced in lfj. Next, we have defined 
F 2n+1 (J; J) = Z- 1 J2 \Jmi\\e- m + e^ E -)[(5{E m - E L )} 2n 

ml 

= l3 2n {[R n Rl + R+Rn]), n = 0, 1, 2, 3, . . . . (20) 

In particular, 

F (J;J) = (J;J), F 1 (J;J) = (JJ + + J + J), F 2 (J; J) = f3([[J + ,n}, J]), 
F 3 (J;J) = f3 2 {[J + ,H][H,J] + [H,J][J + ,n]). (21) 

The functionals ( JIB"]) and ( 120]) are used to generalize all the known inequalities used in the 
Approximating Hamiltonian method. Here we give the final results: 

A. The generalized Harris inequality 

For all integers n = 0, 1, 2, . . . we have 

F 2n {J\ J) < \F2n + i{J\ J) < F 2n {J- J) + ^F 2n+2 {J- J). (22) 

The Brooks Harris inequality ()14p is recovered when n = 0. 



B. The generalized Plechko inequalities 

The following inequalities were proved to hold for all p, q > 1, such that 1/p + 1/q = 1: 
(2Z)- 1 £ | J ml \ 2 \e-' E ' - e-» E -ME m - E,f 



< i( J; J)» {z- £ l^ r"/ ( 2;_^,> l^ - S,IP"«»} 1/8 ■ (23) 

One of the possible choices of p and q here is even integer q = 2k (hence, p = 2k /(2k — 1)) 
which leads to the set of generalized Ginibre inequalities. 



C. The generalized Ginibre inequalities 

These inequalities read {k = 1, 2, 3, . . . ): 

F 2n (J; J) < ^F 2n+1 (J; J) < F 2n (J; J) + i(J; jf k - 1)/2k [F 2k(2n+1) (J; J)] l l 2k . (24) 

At n = the above set reduces to a symmetric version of the inequalities obtained by 
Plechko [id]: 

(J; J) < \(JJ + + J + J) < (J; J) + \(J\ J) (2fc - 1)/2fc /3(i4; R k ) 1/2k , {k = l, 2, 3, . . . )• (25) 



Hence, in the particular case of k = 1 one obtains the Ginibre inequality ([To 



D. The generalized Bogoliubov Jr. - Plechko - Repnikov inequalities 

These inequalities are obtained from (12"3"|) under the choice of odd q = 2k + 1, hence, 
p = (2k + 1) /2k, k = 1,2,3,...: 

^ 2n+1 (J; J) < F 2n (J; J) + ~(J; J) 2fc/(2fc+1) [F 2(2nfc+n+fc)+1 (J; J)] 1 /^ 1 ). (26) 

At n = these reduce to a symmetric version of the set of inequalities obtained by 
Bogoliubov Jr., Plechko and Repnikov [111 ]: 

^(JJ + + J + J) < (J; J) + i(J; J) 2k,{2k+1) {(3 2k (R k Rt + J R+ J R fc )} 1 /( 2fc+1 \ (27) 

The symmetric version of the inequality due to Bogoliubov Jr. (ITT)) follows from here in 
the particular case of k = 1. 



In (7[ we have shown that under sufficiently mild conditions, each of the generalized upper 
bounds has the same form and order of magnitude with respect to the number of particles 
(or volume) for all the quantities derived by commutations of an intensive observable with 
the Hamiltonian of the system. An application of the generalized inequalities to a quantum 
spin model with separable attraction and the Dicke model of superradiance was given too. 

III. FIDELITY SUSCEPTIBILITY AND GIBBS THERMAL STATES 



Over the last decade there have been impressive theoretical advances concerning the 
concepts of entanglement and fidelity from quantum and information theory (l3|,(l4|, and 



their application in condensed matter physics, especially in the theory of critical phenomena 
and phase transitions, for a review see p, [l5| . These two concepts are closely related to each 
other. 

The fidelity (la, ll7| naturally appears in quantum mechanics as the absolute value of 
the overlap (Hilbert-space scalar product) of two quantum states corresponding to different 
values of the control parameters. The corresponding finite-temperature extension, defined 
as a functional of two density matrices, p\ and p 2 , 



1/2 1/2 



%,P2) = ^Vft'V > ( 28 ) 



was introduced by Uhlmann 18[ and called fidelity by Jozsa 16]. 

Being a measure of the similarity between quantum states, both pure or mixed, fidelity 
should decrease abruptly at a critical point, thus locating and characterizing the phase tran- 
sition. Different finite-size scaling behaviors of the fidelity indicate different types of phase 
transitions. The fidelity approach is basically a metric one and has an advantage over the 
traditional Landau- Ginzburg theory, because it avoids possible difficulties in identifying the 
notions of order parameter, symmetry breaking, correlation length and so it is suitable for the 
study of different kinds of topological or Berezinskii-Kosterlitz-Thouless phase transitions. 

The above mentioned decrease in the fidelity J r (pi,p 2 ), when the state p 2 approaches 
a quantum critical state px, is associated with a divergence of the fidelity susceptibility 
Xf{pi) which reflects the singularity of .F(pi,p 2 ) at that point. The fidelity susceptibility 
Xf(pi), which is the main objects of this study, naturally arises as a leading-order term 
in the expansion of the fidelity for two infinitesimally close density matrices p\ and p 2 = 
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Pi + <5p. For simplicity, in this section and thorough the rest of the paper we set in (JI]) 
|A|(zMt + u*A) = hS, where h = v = u* is a real parameter. Following our study [6], we 
consider the one-parameter family of Gibbs states 



p{h) = \Z(h)}- 1 exp[-/3T + (3hS], 



(29) 



defined on the family of Hamiltonians of the form H(h) = T — hS, where the Hermitian 
operators T and S do not commute in the general case, h is a real parameter, and Z(h) = 
Tr exp[— @T+(3hS] is the corresponding partition function. There are two natural definitions 
of the fidelity susceptibility at the point h = 0, depending on the type of approach to that 
point: 

(i) The one-sided fidelity susceptibility at the point h = in the parameter space is 
defined as (see e.g. 



Xf(p(0)) ■= Jim 



-2\nF(p(0),p(h)) _ d 2 F(p(0),p(h)) 



h^O 



h 2 



dh 2 



(30) 



h=0 



Starting from this definition, in our work [6] we derived a spectral representation for xf(p(0)) 
which was used to obtain lower and upper bounds on the fidelity susceptibility. 

(ii) One may consider also a symmetric two-sided definition of the fidelity susceptibility 
at the point h = x G R, given by 

-2 In F(p(x - h/2),p{x + h/2)) d 2 F(p(x - h/2),p(x + h/2)) 



Xf\p(x)) := lim- 

h— >0 



h 2 



dh 2 



h=0 

(31) 



To avoid confusion, we point out that for mixed states the definition of the fidelity 
susceptibility (13"0~|) . based on the Uhlmann fidelity ( 12 8p . differs from the one derived in 



201 ] (see also J5|) by extending the ground-state Green's function representation to nonzero 



temperatures, even both have the same T = limit. This fact has been pointed out in 2l), 



see also our discussion in 



6|. Along with the statistical mechanical notion of susceptibility, 



20 



22], thus 



the quantity (1301) is known also as "Bures metric over the thermal states" 
introducing the geometric approach in the field. 

To proceed with the calculations when the operators T and S do not commute, one has 
to consider a convenient spectral representation. To simplify the problem, we assume that 
the Hermitian operator T has a complete orthonormal set of eigenvectors |n), T\n) = T n \n), 
where n — 1,2, ... , with non-degenerate spectrum {T n }. 
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In the one-sided case (130|) . the following spectral representation was obtained 



ft 



l^|(m|p'(0)|n)| 2 



2^ Pm + Pn 

m,n 

P 2 ^ Pn-pm \(n\S\m)\ 2 1 2 d2 

Here X mn = (3(T m —T n )/2, (• • • ) denotes the Gibbs average value at /i = 0, 55^ = S d —(S d ) , 
where S d is the diagonal part of the operator S, so that 

((^ d ) 2 ) := Y,P m HS\m) 2 - (S) 2 . (33) 

in 

Representation (132]) was the starting point for the derivation of inequalities involving 
macroscopic quantities, like susceptibilities and thermal average values. Note that the first 
term in the right-hand side describes the purely quantum contribution to the fidelity sus- 
ceptibility, which vanishes when the operators T and S commute, while the second term 
represents the "classical" contribution. 

By comparing definitions (130]) and fl3Tl) in the non-commutative case, one sees an essential 
difference: in the one-sided case the zero-field density matrix p(0) is diagonal in the basis 
spanned by the eigenfunctions of the operator T, (m\p(0)\n) = p m (0)(5 mjri , m,n = 1,2, ... , 
while in the two-sided case (J3T1) both p(—h/2) and p(h/2) have non-diagonal elements. 
Hence, it is not obvious that the two-sided fidelity susceptibility will have the same spectral 
representation ( 1321) and, the following from it, lower and upper bounds. This observation has 
motivated us to give an independent derivation of the spectral representation for Xf \P\fy)- 

IV. DERIVATION OF THE SPECTRAL REPRESENTATION FOR THE 
SYMMETRIC TWO-SIDED FIDELITY SUSCEPTIBILITY 



For the family of density matrices of the form (1291) we consider the slightly more general 
case of a two-sided fidelity around a point iGft 

J r (p(x-y),p(x + y)) = Tryjp 1 / 2 (x-y)p(x + y)p 1 / 2 {x-y), y G R. (34) 

Assume that {\n), n — 1, 2, . . . } is a complete set of eigenvectors of the Hamiltonian H(x) = 
T — xS, with eigenvalues E n (x): (T — xS)\n) = E n {x)\n). In this basis the density matrix 
p(x) is diagonal too, which implies (m\p(x)\n) = p n (x)8m,n, for all m, n = 1, 2, . . . . 
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In order to calculate the symmetric two-sided fidelity susceptibility at y = 0, we need the 
0(h 2 ) term in the expansion of (1341) in powers of h — > 0. Consider first the expansions 

p(x ±y) = p{x) ± yp'(a;) + -y 2 p"(x) + r±, (35) 

where r± are bounded operators of the order 0(h 3 ). From the normalization condition for 
the density matrices, Trp(x ± y) = Trp(x) = 1, it follows that 

Trp'(x) = Trp"(x) = 0. (36) 

To calculate the square root of p{x — y) up to the order 0(h 2 ), we set 

p 1/2 (x ±y)= p 1/2 (x) ±A + B + C ± , (37) 

where A, B, and C± are bounded operators of the order of y, y 2 , and y 3 , respectively. Then, 
by comparing the squared of ( 1371) with the expansion ( 1351) for p(x — y), we obtain 

yp'{x) = p 1/2 (x)A + Ap 1/2 (x), 
l -y 2 p"{x) = A 2 + p l '\x)B + Bp x ' 2 {x). (38) 

Hence, we find the matrix elements of the operators A and B: 

H^l»> = y 1/2, vt 1/2 -, v ^ 

Pm [X) + Pn [X) 

, i „, v (m|A 2 |n) 1 o (m|p"(x)|n) . „ . 

(m|S|n) = ~ 1/2 ,' i/2 M + / iVuVn " (40) 

Pm (a?) + Pn (Xj Z p m (Xj + pn \X) 

Next, with the aid of expansions (135]) and (l37j) we evaluate, up to order 0(y 2 ), the product 

p 1/2 (a; - y)p(3 + y)p 1/2 (x - y) = p 2 (x) - Ap 3/2 (x) - p 3/2 (x)A + yp 1/2 (x)p'(x)p 1/2 (x) 

+Ap(x)A - yAp'(x)p 1/2 (x) - yp l/2 (x)p'(x)A + \y 2 p l/2 (x)p\x)p ll2 (x) 

+Bp 3 ' 2 (x) + p 3/2 {x)B + 0(y 3 ). (41) 



Next, following the standard scheme proposed in [23(, we set 



yJpV 2 (x - y)p{x + y)p 1 / 2 (x - y) = p{x) +X + Y + Z, (42) 

where X, Y, and Z are bounded operators of the order of y, y 2 , and y 3 , respectively, and 
compare the terms of the same order of magnitude in the right-hand sides of (14T1) and the 
squared of ( 1421) . Within the order 0(y) we obtain 

p(x)X + Xp(x) = -Ap 3 ' 2 (x) - p 3 / 2 (x)A + yp 1 / 2 (x)p'(x)p 1/2 (x) (43) 
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By taking matrix elements with the eigenvectors of H(x) and using expression (!39j) for 
(m|^4|n), we find 

[PmO) + p n {x)]{m\X\n) = 

-V i/ir'f^ll ; W 2 {*)+PT{*)] +ypll 2 (x)p^(x)( m \p'(x)\n) = 

Prn {X) + Pn {X) 

-y{m\p'{x)\n)[p]l\x)-pH\x)f. (44) 

Therefore, 

H*» = -y ^V^ly ™^ - PT{x)?. (45) 

Pm\X) T" pn\X) 

One important consequence of this equality is the vanishing of the diagonal elements of X, 
hence TrX = 0. 

Next, within the order 0(y 2 ), from (|4"T|) and the square of expression (142]) we obtain 

X 2 + p(x)y + yp(x) = Ap(s)A - yAp(x)p 1/2 (x) - yp 1/2 (x)p'(x)A 
+-y 2 p 1/2 (x)p"(x)p 1/2 (x) + Bp 3 / 2 (x) + p 3/2 (x)B. (46) 

For our needs it suffices to take the diagonal elements of the above equality 

(n\X 2 \n) + 2p n (x){n\Y\n) = (n\Ap(x)A\n) - yp l T l 2 (x)(n\Ap\x) + p'(x)A\n) 
+^y 2 p n (x)(n\p"(x)\n) + 2p 3 J 2 {x)(n\B\n), (47) 

and evaluate Try. Since the calculations are rather involved, we present them in some 
detail. Taking into account Eqs. (|39|) . ( 140]) and ( 145]) . we obtain the expression 

[ P U 2 (x)-pl/ 2 (x)]* 



TrF = ^2(n\Y\n) = --y 2 ^ \( n \p( x )\ m )\ 



2 j^ n ' ' Pn(x)[p m (x) + p n (x)] 

\(n\p'(x)\m)\ 2 p m (x) 2V ^ \(n\p'(x)\m)\ 2 



2 



i 1 ? 2V^ \\ n \P\x)\m)\ p m {x) _ 2ST _ 

^'L. / Nf 1/2/ n , 1/2/ m, y Z^ 1 



2 m,n P„(x)[pm /2 (x) + Pn /2 (x)] 2 m ^ p„ /2 (x) [p^(x) + p« /2 (x)] 

+1 , 2T r/(x) - 1 , 2 V I^Ml^E . (48) 

2 2 ^[p^w+py 2 ^)] 2 l ; 



iii.il 



Now we note that the forth term in the right-hand side of the above equality is zero due to 
(136]) . The sum of the second and fifth terms yields 

1 2 y> \{n\p\x)\m)\ 2 [p m (x) - p n (x)] = 1 2 y \(n\p'{x)\m)\ 2 [pll 2 (x) - pj^jx)} 

2 m,n pn(x)[pU 2 (x)+p 1 n /2 (x)} 2 2 p n (x) [p^ /2 (x) + p^ /2 (x)] 
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Next, by writing [pU 2 - pi/ 2 ] 4 = (p m + Pn) 2 - ±pmPn{p m + Pn) + ^pmpn, we split the first 
sum in the right-hand side of (j4"5j) into a sum of three terms: 



-y^lH^Nf 



2 

m,n 



4:pU 2 (x) 4p m (x) 



(50) 



Pn(x) p 1 J 2 (x)[p m {x) + Pn {x)} \Pm(x) ~ p n (x)f 

By adding up (|49p with the third term in (|48p and the first term coming from ( 15 Op we obtain 

ojy \(n\p'(x)\m)\ 2 

Z!J Z^ 1/2/ U 1/2/ x . 1/2/ si" ^ 0i ^ 

m,n Pn W [P™ (x) + Pn (x)J 

Remarkably, the sum of the above result with the second term coming from (J5"0~l) vanishes: 

2y 2 V IHp^)H| 2 [pm /2 (x)-pn /2 (x)] =Q / 52 x 

m,n [P™ 2 (%) + Pn 2 (x)][p m (x) + p n (x)] 

Therefore, we are left with the contribution of the last term in (1501. hence 



\{n\p'(x)\m)\ 2 p m (x) _ 2V^ \( n \p'( x )\ m )\ 



2 



Try = -2v 2 Y^ iv"ip wi'»/i p™w _ _ 2 y^ ww^ c 53 n 

^^ [p m (x)+p„(x)] 2 ^ ^ [pm(x) + Pn (x)Y l J 

Thus, taking trace of both sides of equality (l4"2"j) . we obtain the two-sided symmetric 
fidelity (jSJ) up to the order 0(y 2 ): 



|(n|p'(x)|m)| 



2 



m,n L v y ' \ /j 

Hence, by setting here y = h/2, from the definition (J5T|) of the fidelity susceptibility, it 
follows that 

At x = this expression reduces exactly to the spectral representation (152]) found in our 
work [6] for the zero-field one-sided fidelity susceptibility. 



V. LOWER AND UPPER BOUNDS ON THE FIDELITY SUSCEPTIBILITY 

Bounds on the fidelity susceptibility follow by applying elementary inequalities for 
(xcothx) -1 to the summand in expression (1521) . In our paper [6j an upper bound on Xf(p) 
was obtained in the transparent form 

Xf( P )<^(SS;8S) , (56) 



13 

where (5S; SS) is the Bogoliubov-Duhamel inner product of the self-adjoint operator 5S 
with itself in the Gibbs ensemble with Hamiltonian H(0) = T. Note that the right-hand 
side of the above inequality is proportional to the initial thermodynamic susceptibility: 

N d 2 f[H(h)} , N 

(5S; 5S) = - — — \ h=0 := —xn, (57) 

where f[H(h)] is the free energy density of the system described by the Hamiltonian H(h) 
and xn is the susceptibility with respect to the field h. 

On the other hand, by applying to the spectral representation for the fidelity susceptibility 
( 152|) the elementary inequality (xcothx) -1 > 1 — (l/3)a; 2 , we have obtained the following 
lower bound 

Xf(p)>^(SS;5S)o-^([[S,T},S}) . (58) 

The quality of the derived upper and lower bounds was tested in the simplest case of 
a single spin in external magnetic field, subject to a transverse-field perturbation. Finally, 
these bounds were applied to two many-body quantum-mechanical models: the single impu- 
rity Kondo model and the Dicke model of superradiance. In conclusion, our lower ( 158]) and 
upper ( 156|) bounds indicate that for the detection of a second order phase transition, with 
diverging in the thermodynamic limit susceptibility, the fidelity susceptibility per particle 
Xf/N is as efficient as the usual susceptibility x- This conclusion is in conformity with 
the commonly accepted view that quantum fluctuations are dominated by the thermal ones 
when T c > 0. However, one should keep in mind that our results were derived under rather 
restrictive conditions on the spectrum of the Hamiltonian. 

VI. CONCLUDING REMARKS 

The infinite sets of generalized statistical mechanical inequalities presented in Section 
II provide upper bounds on the difference between the quadratic fluctuations of intensive 
observables expressed in terms of the corresponding Bogoliubov - Duhamel inner product 
and Gibbs average values of their commutator with the Hamiltonian. Such bounds are 
used, e.g., in the majorization techni que developed by Bogoliubov Jr. for the needs of the 



Approximating Hamiltonian method 10]. A survey of inequalities used to solve problems 
arising in the Approximating Hamiltonian Method, along with their generalizations, is given 
in our paper |7|. The results are illustrated by two types of exactly solvable model systems: 
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one with bounded separable attraction and the other describing interaction of a boson field 
with matter. 

In Section III, some subtle points in the definition of the thermal fidelity susceptibility are 
discussed. The concept of fidelity susceptibility naturally appears as the fidelity's leading 
term in the perturbation expansion with respect to the infinitesimal deviation from a par- 
ticular point of the parameter space. So it is possible to have different definitions depending 
on the way this point is approached. Two definitions are commonly considered: a one-sided 
second derivative with respect of the parameter distinguishing the two density matrices, or 
a symmetric two-sided derivative. In the next Section IV, it is shown that the final result 
for the spectral presentation of the fidelity susceptibility does not depend on which of the 
two definitions is used. While it seems intuitively reasonable, it is nevertheless important 
to prove this statement. 

In Section V, we have presented bounds on the fidelity susceptibility, a notion from the 
information theory, which are expressed in terms of quantities from the statistical mechanics, 
thus emphasizing connection points between these two disciplines. An additional reason that 
stimulates this line of consideration of information-theoretic quantities is that the experi- 
mental setup for measuring thermodynamic quantities is well developed. Thus, estimation 
of metric quantities with the aid of thermodynamic-based experiments seems to be very 
appealing. Note that the fidelity susceptibility reduces down to the Fisher information (for 



details see, e. g., 22|, |24j and references therein) which provides another line of applicability 
of our results to quantum estimation theory. 

We have shown that as far as divergent behavior in the thermodynamic limit is considered, 
the fidelity susceptibility xf and the usual thermodynamic susceptibility x are equivalent 
for a large class of models exhibiting critical behavior. It remains for the future to study the 
effect of the degeneracy of the ground state, especially of a macroscopic one, on the upper 
and lower bounds for the fidelity susceptibility. 
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